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Summary. Abstract We show that the stochastic background of gravitational waves, produced in the 
early cosmological epochs, strictly depends on the assumed theory of gravity. In particular, the specific 
form of the function f(R), where R is the Ricci scalar, is related to the evolution and the production 
mechanism of gravitational waves. Using a neural network algorithm which only requires non-Gaussian 
nature and independence of the input signals we conclude that, in order to detect a CSB signal, the 
interferometric sensitivity of detector like VIRGO will be improved. 
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1 Introduction to the Cosmological Stochastic Background of 
Gravitational Waves (CSB) 

In 1954 Penzias and Wilson, f 1], using an antenna for phone communication, discovered an isotropic 
signal that today is well-known as the Cosmological Microwave Background Radiation (CMBR). Today 
it is commonly accepted as the electromagnetic relic signature of the Big Bang: it has the characteristics 
of a black-body radiation at 3 K. In the same way there would be a relic graviton background originated 
in the Big Bang and the detection of this gravitational signal should provide more information on the 
very early stages of the Universe like the Inflation. 

It is well known that the inflation is an epoch in which the Universe was characterized by an ac- 
celerated expansion, likely due to a phase transition of a scalar field. In particular it was shown by 
Abbott and Wise, |2|, that the perturbation spectra is the signature of some vacuum fluctuations in the 
early Universe and these fluctuations could have produced also gravitational radiation, providing a relic 
signal of Gravitational Waves. The mathematical transposition of this process can be given in terms of 
/(i?)-theories, simply using the formalism of the conformal transformations. In this way, following the 
work of |3|, it would be possible estimate a similar signal in our epoch. Moreover from the amplitude 
of the CSB-GW signal it is possible estimate the average amplitude of the primordial fluctuations, so 
the detection of a relic signal of GW is more than a simple challenge. 



2 Characterization and detection of a CSB of GW signal 

In order to simulate a possible CSB-GW signal, we must account for several characteristics of this 
hypothetic signal. This signal in fact should be isotropic, stationary and unpolarized so that its principal 
property is the frequency spectrum. There are several quantities that describe in a more or less precise 
way this signal: the spectral density S /,(/), the characteristic amplitude h c (F) and the energy density 
per unit of frequency Qowif) = — f^ W t , where / is the frequency of the GW signal, |4|. All of these 
quantities are related with each other, but in order to understand the effects on a detector we need to think 
in terms of amplitude so, using the TT-gauge and making a sort of time averaging of the gravitational 
perturbation detected, h ab (t), we obtain the important results that: 

<h ab (t)h ab (t)>=2 f d(\ogf)h]{f). (1) 

Jf=0 

We can also obtain the contribute in terms of energy density of the gravitational perturbation. If we 
make a spatial averaging over several wavelenghts we obtain the density of stochastic GWs paw, which 
relates directly with the characteristic amplitude. Now dividing by the critical density of the Universe, 
we obtain: ^ 

Oowif) = ^fh]{f\ (2) 

which represents a numerical contribution to the total density parameter of the universe Q. Note that, 
as we will see, h c is very small, so this contribution is negligible in the well-known cosmological 
consideration about the total density of the Universe. 

Actually, h c (j) is not yet the most useful dimensionless quantity to use for the comparison with 
experiments. In fact, any experiment involves some form of binning over the frequency. In a total ob- 
servation time T , the resolution in frequency is Af = l/T , so it is convenient to define: 

W^/) = W)fy) . (3) 

and using l/(lyr) ^ 3.17 x 10 s Hz as a reference value for hoQcwif), we finally have the expression 
for the characteristic amplitude: 
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h c{ f,Af) - 2.249 x 10- ffif (*2^f f *L f . (4 ) 
u 7 \ f ) \ 1CT 6 j ^3.17x lQr»Hz) 

This quantity matches directly with several constraints given by some experiment by satellite as 
COBE and WMAP, but also from the Big Bang Nucleosynthesis theory, for which we would have 

K(f)< 2.82 x 10- 21 |^J. (5) 

Moreover another constraint could be given by the Sachs-Wolfe effect for which 




where ST corresponds to the anisotropies of the temperature in the CMBR. 

Now the detected signal is directly connected with the gravitational perturbation by s(t) = iy b h ah (i), 
where D ab is the detector tensor which depends on the directions of the GW-interferometer arms. How- 
ever the output signal will present an additional component due to the interferometer noise so our sim- 
ulated signal will be S (t) = n(t) + s(t), where n(t) is the noise signal component. Using a pre-compiled 
MATLAB algorithm, |5], we reconstructed n(t) for a GW-interferometer like VIRGO, while for s(t) 
we considered several oscillated (sinusoidal-like) signal; this because we are interested in the simple 
detection, not in the rivelation of the wave form, of the CSB signal. 



3 Signal Analysis 

With this in mind we recognize that the reconstructed signal is a chaotic signal, so from the observation 
of a dynamical variable we want to find a method to reconstruct the phase space, or the number of 
variables that compose the starting signal. 

For this reason we used the embedding theorem, (6j, which states that, given an N-point time series, 
X\ , x 2 , Xjy, the phase space vectors are reconstructed in order to get a signal matrix as 

X[ Xl ••• XN-(m-\) T 

X\+t X2+r ' ' ' -^V-(m-l)r+T /j\ 

^Xl+( m -l) T X2+(„,-l) T ■ ■ ■ Xjv 

where r, the time delay, and in, the minimum embedding dimension, or the dimension of the recon- 
structed phase space, are two unknown parameters to determine. 

In order to determine t we minimized the average mutual information 1 T : 

It = £ J P(Xi,x i+T )log2 

that represents the average amount of information about x j+T conditioned by observations x t . For m we 
used a method found by Cao, |7|, which, without going into details, allows to determine the minimum 
embedding dimension. 

The knowledge of t and m allows us to use a statistical algorithm, the Independent Component 
Analysis (ICA), |8], derived from neural networks, in order to going back to the single sources which 
compose the starting signal. If we observe a signal y = Ax, with A mixing matrix and x the unknown 
sources, that satisfies the following conditions: a) they must be independent, b) they must be non- 
Gaussian, then we can construct a combination of the observed signals as: 



P(x,)P(x, + 



(8) 
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s = w T y = w T Ax = z T y. (9) 

Now for the Central Limit Theorem s becomes less Gaussian when equals any one of the y. In general s 
is more Gaussian of any y, so we need an estimator of non-Gaussianity. This is given by the Negentropy, 
which states for negative statistical entropy, J(s) oc [EG(s) - EG(v)], where E is the expectation value 
of the quantity in consideration and v is a null-mean and unitary-variance Gaussian variable. In this 
way we require a sort of maximization of the Negentropy and we obtain it using an iteration algorithm, 
named FastlCA. 

The algorithm FastlCA consists of 2 steps. The first one is an iteration cycle for the determination 
of a single unit, i.e. a single vector component of the inverse matrix w. The cycle is the following: 

• random choice of the starting vector w 

• w = E{yg(w T y)j -E|g'(w T y)|w 

• w = w T / || w T || 

• if there isn't a convergence, do repeat the cycle 

where convergence stands for new and old values of w that point in the same direction. The second 
step is the determination of the remaining components, for which we used a two-steps Gram-Schmidt 
algorithm: 

p 

Wp+l = ^p+l - ^ W^WjWj, (10) 

;=i 

and 

W p+ i = W pt i/Jw^W^i, (11) 

that therefore allows us to recover w and therefore the mixing matrix A, from which, and from y, we 
come back to the single components x of our starting signal S(t). 

So we have constructed an algorithm for the estimate of the independent component of an hypothet- 
ical CSB-GW VIRGO-detected signal. In this way, using our procedure, we can determine a minimum 
detectable threshold for the CSB-GW signal, giving to the gravitational signal s(t) an oscillatory be- 
haviour varying with the frequency. The results of our analysis are shown in figure [3] and in the table 
[T] We immediately see that the maximum sensitivity is obtained for signals of frequency in the range 
(100- 1000Hz), but more important we note that our algorithm works very well with this type of signals 
also if the identification of the estimated signals were done with a graphical recognition. 



Table 1 Results of the simulations 



Frequency (v) Amplitude h(t) 



0.3 


io- 19 


0.7 


io- 19 


1.3 


5 x 10~ 20 


2.5 


5 x 10- 20 


6.5 


5 x 10- 20 


100 


5 x 10- 20 


200 


IO" 20 ~ 5 x 10~ 20 


400 


io- 20 


650 


5 x 10- 20 


1500 


5 x 10- 20 ~ 10- 19 


3000 


io- 19 
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Fig. 1 Minimum detectable threshold by an interferometer VIRGO-like for an hypothetical CSB signal 
of GW, with frequency increasing. The maximum sensitivity is obtained for signals of frequency in the 
range (100 - 1000#z). 

4 Conclusions 

From our analysis we conclude that it is necessary to constraining the theory for a Cosmological 
Stochastic Background of Gravitational Waves and we advise to deepen the relations among the Sachs- 
Wolfe effect and the gravitational perturbations in f(R) theories. This could be a starting point for the 
theoretical analysis of the detection with interferometer of a gravitational signal coming from the relic 
background originated from the Big Bang. 

Moreover also the application to VIRGO-detected signals of the ICA methods needs to be im- 
proved, in particular using a time series longer, a number of the sources more than 2 and possibly the 
utilization of a less-human method for the identification of gravitational signals. The last, but not the 
least, we suggest a frequency analysis in order to use directly the energy density for the study of the 
connections with the Sachs-Wolfe effect. 
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